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Abstract 

We examine the derivation of the spin-orbit force from the Skyrme model. 
We find substantial agreement with phenomenological potentials. We show that 
a systematic and simultaneous treatment of all components of the kinetic energy 
introduces a minus sign relative to previous calculations for the iso-spin inde- 
pendent part of the interaction. We also explain the smallness of the iso-spin 
dependent part of the spin-orbit potential. 

The challenge of obtaining the A^A^ interaction from a QCD inspired starting point 
has recently been given new impetus by results in the Skyrme model [0, 0. Careful 
treatment of the non-linear nature of the model and of quantum corrections leads to 
significant central attraction, in substantial agreement with phenomenology. There is 
even a strong short range repulsion, though its detailed description certainly lies outside 
the scope of the Skyrme model with pions only. Combined with the well known Skyrme 
treatment of one pion exchange and the tensor force, all this yields a satisfactory picture 
of the static part of the NN interaction that gives a good description of the phase shifts. 
The next step is to consider the non-static parts of the interaction. Here we consider 
the spin-orbit force. The dominant part of the phenomenological spin-orbit interaction 
is iso-spin independent. Attempts to reproduce it from the Skyrme model have given 
the correct order of magnitude, but the wrong sign |P, In this note we point out 
that a careful treatment of the conversion from velocities to canonical momenta in fact 
reverses that sign and hence yields a spin-orbit interaction in substantial agreement 
with the data. We also show that a small iso-spin dependent spin-orbit force emerges 
naturally. 

The principal ingredients in the successful static A^A^ interaction calculations are 
numerically "exact" solutions for the baryon number two {B = 2) system at finite 
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separation, and quantum corrections implemented through the Born-Oppenheimer ap- 
proximation. The exact results differ considerably from earlier ones based on the 
simple product Ansatz largely because the static interaction energy arises as the rela- 
tively small difference of two large numbers. Cancellations must therefore be treated 
with care. Although a full treatment of the non-static terms should also be done ex- 
actly, there are no corresponding cancellations in this case, and one can hope to get 
insight from the product Ansatz. It is for that reason that the wrong sign for the 
spin-orbit interaction obtained in earlier product Ansatz calculations was particularly 
disappointing ^. 

Since there has not been a complete study of the kinetic energy of the Skyrme 
model, even though pieces can be found in the literature P, ||, we decided to 
undertake such a study. A part of our purpose was to understand the puzzle of the 
spin-orbit sign. In this note we sketch the steps in our evaluation, and show that there 
is a simple solution to the sign problem. 

The Skyrme model is a non-linear field theory that can be realized in terms 
of an SU (2)- valued matrix field U, with Lagrangian density 

C = ^TT[d^U{x)d>'U^{x)] + ^Ti[U^d^U, U^d.Uf + ^[Tr(f/ + f/t) - 4]. (1) 

The model is covariant, as well as invariant under global SU(2)-rotations that are 
identified with the iso-spin symmetry. The only parameters in the model are the 
pion mass, m^, the pion decay constant, /^r, and the dimensionless parameter, g. 
It is customary not to take /tt from experiment but rather to adjust /^r and g to 
give the nucleon mass and the nucleon-delta splitting correctly. We take f-,^ = 6.45 
MeV, g = 4.48 and use = 138 MeV. It should be kept in mind that all the 
results on the static A^A^ force and on the spin-orbit interaction are based on only 
these three parameters. As was discovered by Skyrme the model has a topologically 
conserved quantum number, which is identified as the baryon number, B. The U 
field is interpreted as a combination of a scalar a field and an iso-vector pion field, 
U = {a + ir ■ Ti)/ fn- The a field is not an independent physical field due to the 
unitarity constraint on U. 

The standard time-independent solution to the classical field equations for B = 1 
is the defensive hedgehog, where the pion field points radially outward, 

f/i (r) = exp(ir ■ f/(r)). (2) 

The baryon number of this state is given hj B = (/(O) — /(oo))/7r = 1. This solution 
breaks translational invariance, as well as the 0(4) spin-iso-spin symmetry, but the 
sum of spin and iso-spin is still conserved. If we perform a global SU(2) iso-rotation 
on the state, 

U,if\A) = A^U,ir)A, (3) 
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we obtain a state of the same energy. 

For the B = 2 system, we will use the product Ansatz. This Ansatz makes use of 
the fact that the product of two B = 1 solutions has baryon number two. The most 
general Ansatz we can construct from two hedgehogs consists of the product of two 
separated and rotated hedgehogs, 

U2ir\RAB) = A^Ui{f-R/2)AB^Ui{r + R/2)B 
= U2{r\RCD) 

= D^C^'^^Ui{r-Rl2)CUi{r + R/2)C^'^^D. (4) 

In the last line of (1) we have introduced the matrix that describes the rigid iso- 
rotation of the whole system, as well as a relative iso-rotation C . When R is very 
large changing C ot D does not change the energy of the solution. For smaller R, 
D still generates a zero- mode (corresponding to broken iso-spin symmetry), but the 
energy will depend on C. The energy is also invariant under spatial rotation, due to 
the conservation of angular momentum J = L + S . 

In order to evaluate the kinetic energy, we make all the collective parameters 
{R, A, B) time dependent, and substitute the product Ansatz (4) into the Lagrangian. 
This leads to an effective Lagrangian containing a kinetic energy that is second order 
in the time derivatives due to the special nature of the Skyrme Lagrangian. Note that 
we have already used Galilean invariance to decouple the center-of-mass motion from 
the relative motion. 

We have evaluated all terms in the kinetic energy. The terms can be expanded in 
tensors in C and powers of R (as well as the velocities, of course). The tensors form a 
basis for the irreducible representations of the 0(4) formed from the four components 
of C. We find only symmetric representations of 0(4) appearing, (cr, 0), with a even. 
Only (7 = (the unit operator), cr = 2 and o" = 4 occur in the expansion. Each of these 
terms is multiplied by a coefficient function that depends only on the size of R. We 
find typically that these functions are numerically largest if cr = 0, smaller if cr = 2, 
and smallest if cr = 4. Note that the cr = 4 operators have zero matrix element in 
two nucleon states. We also find several non-zero coefficients that are not even under 
parity or the interchange of the two Skyrmions - a well known disease of the product 
Ansatz. In the remainder of the discussion we ignore these terms. The full expansion 
of the kinetic energy is quite complicated and will be reported on separately Here 
we concentrate on the spin-orbit coupling and only consider terms that are dominant 
in that expansion. 

In order to quantize the kinetic energy, we must first invert the mass matrix to 
obtain the Hamiltonian In other words we perform the standard transformation 
from a Lagrangian, 

L = \Y,q,M,,q^~V{q), (5) 
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to a Hamiltonian, 

H=lY.MM-%p, + Viq}, (6) 

ij 

with Pi = J2j MijQj conjugate to g^. We separate the mass matrix in terms of its 
dependence on the 0(4) label, a, 

M = Mo + M2 + M4. (7) 

As mentioned above these terms decrease rapidly in size with increasing a. We take 
advantage of this fact to drop M4 and write 

M^"^) ^ Mo~^ - Mo'^MaMo"^ + small terms (8) 

The LS coupling already appears in Mq, as an off-diagonal terms coupling the velocity 

V = R to the sum of the two rotational velocities for the individual Skyrmions. If 
we denote the sum and difference of the rotational velocities by u^, we find that the 
dominant contribution of Mq to the kinetic energy is given by 

K = l{2d\{Ry + d^(i?)(cu+' + u'^) + 2ul^{R){v x R ■ cD+)}. (9) 

Here the quantities ,and are functions of R, calculated from the Skyrmion 
configurations. We have that d\ ^^-^ M/4, d'^ 2 A and u^"^ ^zl^ g (M is the mass 
of a single Skyrmion, and A its moment of inertia.) 

The approach in P, H was to replace v by p/(M/2), and a;"*" by S/{2A) to find 

V s rim 

However this procedure is not equivalent to inverting the mass matrix. The matrix to 
be inverted here is basically two-by-two with the spin-orbit coupling coming from the 
off diagonal part. As is well know the off diagonal terms in such a case change sign 
under inversion. If we carry out the inversion from (5) to (6) we indeed find for the 
spin orbit force 

2dldf-{uf'y MA ^ ' 

Thus we find the correct sign for the spin-orbit interaction. In Fig. 1 we show our 
calculation for the iso-spin independent part of the spin-orbit interaction. We show 
two calculated curves compared with the BonnR potential |10[. One calculation uses 



the asymptotic values (large R) for d^ and df for all R and neglects the ^3"^ in the 
denominator (this is close to the spirit of 0]), while the other keeps all the terms in Eq. 
(11) with their full R dependence. Both curves agree fairly well with phenomenology 
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Figure 1: The isospin independent spin-orbit potential as a function of nucleon sepa- 
ration R. The sohd hne is calculated using the first expression in Eq. (11), whereas 
the dashed line represents the result from the second expression. The dash-dotted line 
is the corresponding term in the BonnR potential flO . 
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Figure 2: The isospin dependent spin-orbit potential as a function of nucleon separa- 
tion R. The solid line is our result, and the dash-dotted line the corresponding term 
in the BonnR potential. 



beyond 1 fm, where the product ansatz makes sense, but differ at small distances. The 
R dependent curve is probably the better indicator of what a full calculation will yield. 

The iso-spin dependent part of the spin-orbit interaction comes from the M2 term 
in Eq. (8). It is shown in Fig. 2 and compared with the BonnR potential ||T0|- As 
discussed above it is smaller than the iso-spin independent part, as is the data. It is 
also complicated with many terms contributing. As seen in the figure, it is positive at 
large i?, and oscillates at smaller R. This oscillation comes from a cancellation among 
many terms and is probably quite sensitive to model assumptions. This sensitivity 
to cancellations recalls the static interaction in the central channel, suggesting that 
channel coupling may again be important 

Riska and Schwesinger have added a 6th order term to the Skyrme model, and 
claim to find the right sign for the spin-orbit interaction. The results of this note shows 
that their discussion is faulty. It is not necessary to include a 6th order term in the 
Skyrme model to find agreement between the model and phenomenology for spin-orbit 
coupling. Furthermore their 6th order term calculated as in (6) and (7), will lead to 
the wrong sign of the spin orbit interaction. 

In summary we have shown that the Skyrme model yields the correct sign and order 
of magnitude for the NN spin-orbit interaction. Coupled with our earlier demonstra- 
tion that the Skyrme model correctly accounts for the main features of the static 
A^A^ interaction, one pion exchange, mid-range attraction, tensor force, short range 
repulsion, the result concludes successfully the first steps in the derivation of the NN 
interaction from a starting point based in non-perturbative QCD. 
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